Local random quantum circuits are approximate polynomial-designs - numerical results 
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We numerically investigate the statement that local random quantum circuits acting on n qubits 
composed of polynomially many nearest neighbour two-qubit gates form an approximate unitary 
poly(n)-design [F.G.S.L. Brandao et al., arXiv:1208.0692|. Using a group theory formalism, spectral 
gaps that give a ratio of convergence to a given i-design are evaluated for a different number of 
qubits n (up to 20) and degrees t (t = 2, 3, 4 and 5), improving previously known results for n = 2 
in the case of t = 2 and 3. Their values lead to a conclusion that the previously used lower bound 
that bounds spectral gaps values may give very little information about the real situation and in 
most cases, only tells that a gap is closed. We compare our results to the another lower bounding 
technique, again showing that its results may not be tight. 
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I. INTRODUCTION 

Random unitary matrices have established their place 
as a useful and powerful tool in theory of quantum in- 
formation and computation. From various examples 
of their applications, let us report here that random 
unitaries are used in the encoding protocols for send- 
ing information down a quantum channel [lj, approxi- 
mate encryption of quantum information |2], quantum 
datahiding [2], information locking [2], process tomog- 
raphy J3J, state distinguishability [4], and equilibra- 
tion of quantum states J5HH or some other problems 
in foundation of statistical mechanics [9|. 

But their is a problem with them - they are not very 
favorable from a computational point of view. Why? 
The answer is the following - to implement a random 
Haar unitary one needs an exponential number of two- 
qubit gates and random bits (in other words, to sam- 
ple from the Haar measure with error £, one needs 
exp(A"log(^)) different unitaries). To omit this prob- 
lem, one can construct, so-called, approximate random 
unitaries or pseudo-random unitaries. Using them, an 
efficient implementation is possible. 

An approximate unitary £-design is a distribution of 
unitaries which mimic properties of the Haar measure 
for polynomials of degree up to t (in the entries of the 
unitaries) 1 10-23 1. Approximate designs have a number 
of interesting applications in quantum information the- 
ory replacing the use of truly random unitaries (see e.g. 
ED Ql El E2H21 What is more, other particular 
constructions of approximate unitary £-designs together 
with some applications in quantum physics, have been 
formulated, let us mention here, for example, a recent 
construction of diagonal unitary f-designs If26ll27l . 

In 2009, Harrow and Low [15J stated a conjecture that 
polynomial sized random quantum circuits acting on 
n qubits form an approximate unitary poly(n)-design. 



To give an example supporting their statement, the au- 
thors, also in 2009, presented efficient constructions 
of quantum f-designs, using a polynomial number of 
quantum gates and random bits, for t = O(log(n)) ||23] . 

However, it took some time to verify the conjecture 
from 1 15 1 . At that time, it was already known that there 
exist efficient approximate unitary 2-designs in V(d n ), 
where efficient means that unitaries are created by poly- 
nomial (in n) number of two-qubit gates and the distru- 
bution of unitares can also be sampled in polynomial 
time (in other words that random circuits are approx- 
imate unitary 2-design) ITT3HT81 |28l . In 2010, Brandao 
and Horodecki |21J, made one step further proving that 
polynomial random quantum circuits are approximate 
unitary 3-designs. Quite recently, a break-through has 
been made for the above problem. In 1291 , authors 
has proved that local random quantum circuits acting 
on n qubits composed of polynomially many nearest 
neighbor two-qubit gates form an approximate unitary 
poly(n)-design, settling the conjecture from |15| to be 
affirmative. Their proof is based on techniques from 
many-body physics, representation theory and combi- 
natorics. In particular, one of the tool used to obtain the 
main result was estimation of the spectral gap of the 
frustration-free local quantum Hamiltonian that can be 
used to study the problem (instead of a quantum cir- 
cuit). 

In this paper, we analyze two aspects of their state- 
ment: First, we numerically verify and investigate it, 
calculating spectral gaps for the increasing number of 
qubits (in circuit) n (up to n — 20, in the best case) 
and degree t (f = 2,3,4,5). Previously, exact values 
were known only for n = 2 and t = 2,3. This may be 
of independent interest since, in many-body physics, 
the knowledge of spectral gaps of local Hamiltonians 
is useful in studying many-body systems (see, for ex- 
ample, I5l l30rl32| ). In addition, in the authors ob- 
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tained that for n — 2, gap values for f = 2 and 3 are 
the same; we obtain the matching of all calculated spec- 
tral gap values in the case of t = 2 and 3, which is a 
unique feature for this exact degrees f. Second, in Il29l , 
in order to proof the main statement, a lower bound for 
spectral gaps was derived (independent from the num- 
ber of qubits). However, based on our results, it can 
be concluded that lower bounding of spectral gaps may 
not be tight. We show that there could be a large differ- 
ence in actual values of spectral gaps and correspond- 
ing lower bounds. What is more, our calculations for in- 
creasing values of n and f lead to a conclusion that first, 
the lower bound from [29 j is hard to estimate, second, 
it gives only the limited information about the actual 
value of spectral gaps (that a gap is closed). Simultane- 
ously, we also compare our results to the another lower 
bound, derived in [33 J, and show that predictions ob- 
tained according to it can be better than that from Il29l . 
Nevertheless, also in this case, predictions according to 
that bound can be inconclusive in some cases, showing 
that, in principle, obtaining a good bound is a demand- 
ing task. 

Our paper in structured as follows. In Section |n] we 
start from introducing local (random) quantum circuits 
which, using the formalism of superoperators, we con- 
nect with approximate unitary f-designs. In Section III 
we recall the statement that local quantum circuits ofa 
given length form approximate unitary f-designs. We 
also recall that to verify and check this connection, cal- 
culations of spectral gaps and second largest eigenval- 
ues of local Hamiltonians connected to the problem are 
necessary. Section IV is the key section of this paper. 
We first show, how to connect our problem of check- 
ing the connection between local quantum circuits and 
f-designs (calculating the spectral gaps) with symmet- 
ric groups S(f), where f corresponds to a degree in f- 
design. Then, we present our numerical calculations 
for spectral gaps for different number of qubits n and 
different degrees f which we later compare to results 
obtained using the techniques for lower bounding the 
spectral gaps. 



II. RANDOM UNITARY CIRCUITS AND 
APPROXIMATE f-DESIGNS 

In this section, we present the formalism of local 
random quantum circuits and approximate unitary de- 
are 



signs. We want to point here, that Sec. Inland III 
mainly based on 12T1|29] , so for the full analysis (proofs, 
etc.), we refer to these papers. 

We consider n qubits (so, from now d = 2), and apply 
k steps of a random circuit (random walks on U(d n ). 

Definition 1. (Local quantum circuit) In each step of the 
walk an index i is chosen uniformly at random from [n] and 
then a unitary Uy+i drawn from the Haar measure on V(d 2 ) 
is applied to the two neighboring audits i and i + 1. 



Now, we will briefly introduce the formalism of 
norms that will be used to compare two superopera- 
tors. For X 6 B(C d ) we define the p— Schatten norms 

i 

||X||p := tr(|X|P)f. then for a superoperator A : 

B(C d ) -> B(C d ) we define the p ->■ q induced Schat- 
ten norm as 



|A(X)| 



sup 



|A(X)|| 
XL 



(H!) 



Finally, the diamond norm is defined using the 1 ->1 
norm as 



A 



sup 

d 



) A| 



1H>1- 



(n.2) 



There are several different definitions of 
£— approximate unitary f— designs |21] from which let 
us mention the following. 

Definition 2. (Approximate unitary t-design) Let {u, U} 

be the distribution of unitaries after one step of the random 
walk. Define 



V{d) 



U 0t p(U f ) 0t ii(dU) (113) 



and 



Q m {p) = / U^piU^yuidU), (II.4) 

where uh is the Haar measure. Then k steps of random cir- 
cuit constitutes a e-approximate unitary t-design if 



(II.5) 



Having establish these, we can ask the following 
questions 

• How well a random circuit can mimic a random 
global unitary? (i.e. how fast the related random 
walk converges to the Haar measure) 

• How well a random circuit can mimic 'twirling' 
/ dUU <g> IT <g> . . . ® ® tf + ® . . . <8> U f ? (i.e. 
how fast the random walk converges to so-called 
f-design) 

• More generally: How fast a random circuit con- 
verges to f-designs? 

Recently, in a series of papers, an effort has been made 
to answer them ET1I291 . 



III. LOCAL RANDOM QUANTUM CIRCUITS ARE 
APPROXIMATE POLYNOMIAL-DESIGNS 



In this section we will focus on local random circuits 
and show the extension of their formalism to uniform 
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random circuits. First, let us explore some connections 
between superoperators and operators. For a superop- 
erator G given by 



G(X)=£A k XBl, 

k 



(m.i) 



where t denotes the Hermitian conjugate; we define the 
operator 



(III.2) 



with ★ being the complex conjugate. 

Let X be a normalized operator, such that Tr(XX + ) = 
1. What is more, assume that Q(X) = AX, for a com- 
plex eigenvalue A, i.e., X is an eigenoperator of Q with 
eigenvalue A. Then defining |X) = X <8>I|#), with 



ltf>)=£|fc)<s>|fc), 

k 



(1113) 



it holds that G|X) = A|X), i.e., |X) is an eigenvector of 
G with eigenvalue A. 

A direct implication of this correspondence is that 



12^2 



(m.4) 



where the norm 2 — >■ 2 is the norm, defined as in Eq. 

Now, let us connect the problem of local random 
quantum circuits with many-body physics techniques 
and define a local quantum Hamiltonian (a local Hamil- 
tonian on n D-dimensional subsystems is a Hermitian 
matrix H, acting on (C D )® n ) as 



k 



(III.5) 



where each H k acts non-trivially only on a constant 
number of systems. In [21 , 29 1, the authors, introduce a 
specific local Hamiltonian, with D = d 2t 



n-l 

i=l 



(111.6) 



with local terms — I — Ptf+i/ where I is the iden- 

tity and Pjj+i is defined the following: 



p i,i+l 



mti®("y+i)®VH(du). (m.7) 



For the structure and properties of H„ f, we refer to 
Appendix VIII This Hamiltonian allows to formulate a 
very useful (and crucial for later calculations) theorem 

Theorem 3. Let {fi(dU), U} be the distribution ofunitaries 
after one step of the random walk according to the local ran- 
dom circuit model. Then, using Eq. flll.2fr , define 



G t i* k ,t 



V(d) 



U 



(III.8) 



where 



so fi* k is the k-fold convolution of u with itself. Then, we 
have that 



*k 



J 6 Ul ... Uk ii(dUi)...ji(dU k ) r (111.9) 



J }lH,t\ 



< 1 



A(H„, 



(111.10) 



with A(H n j) being the spectral gap of the local Hamiltonian 
acting on n qubits and jih, the Haar measure on U(2"). 

After a successful estimation of the spectral gap from 
Eq. I pTI.lO) , one can show that 

Theorem 4. [29] Local random circuits of size 
(log C^j log(d)log(t))t 5 n 2 form an e-approximate unitary 
t-design. 



At the end, let us note that Eq. | |III.10 \, can be rewrit- 
ten in a such form that it gives relation between the 
spectral gap and the second largest eigenvalue of an 
operator X 



Fact 5. We have 



|G ,k it - G 



}l H ,t \\co 



(III.ll) 



where A2 is the second largest eigenvalue of G^t. Moreover 
the largest eigenvalue X\ of Gu,t is equal to 1, and the corre- 
sponding eigenprojector is equal to Gu H ,t- 

Thus the problem reduces to analysis of spectral gap 
of the operator G^f Now, it is important to ask the 
following: 

• How a spectral gap depends on number of qubits 
n and the degree of design f? 

In the next section, an answer to this question is pro- 
vided. 



IV. SPECTRAL GAPS - NUMERICAL RESULTS 

In this section, we present our numerical results for 
spectral gaps for different degrees t in f-designs (for 
simplicity, we gather all results in Table |IJ. What is 
more, we (where it was possible) compare our results 
with two lower bounds for spectral gaps: 1) The "lo- 
cal" lower bound obtained in [33 1; 2) The "global" lower 
bound derived in [29 1. To evaluate the spectral gaps, 
the Matlab software and a C + + code have been used. 



A. Local Hamiltonians as a tool for calculating spectral 
gaps 

We already showed that calculations of spectral gaps 
can be connected with the second largest eigenvalue of 
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Gu,t- Now, we will show how to connect calculations of 
spectral gaps (equivalently, second largest eigenvalues) 
with symmetric groups S[t), where t plays a role of the 
degree of £-design; usin g te chniques for local Hamilto- 
nians introduced in Sec. TTlTl 

At the beginning, let us remind that our Hamilto- 
nian is of the form H n j = YIlZx h-ij+V with local terms 
hi = I — Pij+i, and the notation as in Sec. Ill 

Let us consider sup erope rators associated with pro- 
jectors Py+i from Eq. ( |IIL7) 



U(d) 



U^X(Ul l+i r^u(dU) (IV.l) 



Now, we can find, as a consequence of the Schur- 
Weyl duality |34] that all operators X invariant under 
action of V{ i+\ can be written as a sum of permutation 
operators V n n i+ iy permuting t copies of the Hilbert 

space Hfj +V Because V n m + ij = V n/i (g V n j + i it follows 



that P iri+t c Pi 
the expression 

Pi = 



P,+l, where operators P, are given by 



U 



U(d) 



(IV.2) 



Thanks to the above consideration we can deduce that 
operator X can be identify with the operator Gu,t from 
Fact [5] and written according to the formula 



G }i ,t = P\,2 <8> h ® . . • <E> In + h <8 Pz,3 «l 14 ' 
+ . . . + li <g> I 2 <g> . . . <g> P n _i „. 



(IV.3) 



In a subspace spanned by permutation operators acting 
on Hilbert space Ti® we are able to construct operator 
basis which is orthogonal in the Hilbert-Schmidt scalar 
product (see Appendix |VI[|. Using this basis we can cal- 
culate two biggest eigenvalues of operators X which are 
necessary to know the spectral gap. 



B. Lower bounding a spectral gap 

Herewith, we present two methods for lower bound- 
ing the spectral gaps that we will use later to compare 
to values of spectral gaps obtained numerically. 

1. "Global" bound: 

Lemma 6. 129] For every integers n (number of qubits) and 
t (t in t-design), with n > [10Zog(f)], the spectra] gap 
A(Jf„ ; t), can be fower bounded as follows: 



A(H„, t ) > 



A(H 



\2(log(d))-Uog(t)-],t 

8(log(d))-Vog(t) 



(IV.4) 



with d being the dimension of the Hilbert space and \a~\ the 
cetiing function, so the function which gives the smallest in- 
teger > a. W 



2. "Local" bound: 

Lemma 7. j|33"P For every integers n and t, the spectral gap 
A(H Ur t) of local Hamiltonian can be lower bounded in the 
following way: 



A(H M ) > 



*A(ff M ) 



(TV5) 



where A(Hj. t ) is the Hamiltonian restricted to k qubits: 



Hi- 



Let us explain now, why we used the terms "global" 
and "local" to describe these two bounds. From Eq. 
I TVA\ it can be noticed that to lower bound a given spec- 
tral gap, one needs to compute only one spectral gap, 
the one that is given by t and the dimension of the sys- 
tem (but we consider qubits only, so the value d = 2 
is set). The bound is global in this sense that it does 
not change with the number of qubits - it remains con- 
stant for an arbitrary length of a quantum circuit, so it 
is independent from qubits number. On the contrary, 
the second bound (from Eq. ( |IV5) ) has totally reverse 
property. There, the gap for a higher number of qubits 
implies better precision of the bound. 

Of course, as we will see later, both bound have pros 
and cons. For example, the "global" lower usually lower 
bounds the spectral gap in a harsh way (only giving an 
info that a gap is closed) and, what is more, one can 
observed that for big values of t, to estimate the lower 
bound, one needs to know a value of the spectral gap 
for a big number of qubits to lower bound the absolute 
value of the gap. Keeping in mind that calculating spec- 
tral gaps is, in principle, computationally hard problem, 
the effectiveness of the "global" bound is limited. The 
good thing about it is that it is always positive (in the 
sense of value). On the other hand, it is not true for the 
second bound. For small numbers of qubits, the "local" 
bound could be negative and it means that one needs 
to increase its precision by calculating the lower bound 
for a gap for a bigger number of qubits. What is more, 
the value of bound sometimes fluctuates - so the value 
of the bound for a number of qubits k could possibly be 
worse than that for k — 1 qubits. The advantage of this 
bound is the following, when it is positive, its value is 
usually closer to the exact value than this predicted by 
the bound from Eq. ( |IV4) . 



C. Method of Calculations 

In this section we will present the reader some meth- 
ods of calculations used in this paper. Notation is 
mostly taken from [21 j. 

We know that for an arbitrary subspace of operator 
space we are able to construct basis of operators which 
is orthonormal in the Hilbert-Schmidt scalar product. 
For this construction we use linear combination of 
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nonorthogonal permutation operators acting on H®*: 

Rk= E hn(t)V n . (IV.6) 

7teS(t) 

Using ope rators K/t we can rewrite operators Py+i 
from I III. 7 1 and Pj form ( IV. 2 1 in a form 



(0\/p(0i 

k I 



To illustrate the action of operators R k 
t x n lattice (see Figure 



(v+i) 



(IV.7) 
consider a 



systems from z* and (i + l) th column. Finally to obtain 



1\, then _r[' /I+1 ^ acts jointly on 



th 



: 



i i+1 



11 



FIG. 1: (Color online) The figure presents a t x n lattice. 
We have n columns and in every column there are t qubits 
on which permutations from S(f) act. In a blue rectangle, 
we marked 2 th and (i + l) th column on which the operator 
+1) acts. 



a representation of an operator X from equation [ IV.3[ 
in our product basis i.e. 



(l) 



kl 



< B) ><< } 



R^\, (IV.8) 



where k — (k\ . . .k n ) and 1 = (l\ . . . l n ) are multiindices, 
we have to express operators P,- (+ i in terms of a product 
basis. For this purpose we use argumentation from [21], 
having: 



Rl' 2 = E r suRs®R 



(TV.9) 



(k) 

where r su are some coefficients which we want to 

know. In this paper, to calculate numbers fs« we use 
the Schur basis (see, for example, [35]). Then every 

1 2 

operator R,' corresponds to a linear combination of 

Ejy <8> Emn/ where operators Ef, form an operator basis 

in a given invariant subspace of H®* labeled by a. Now 
we see that an index k in Eq. | |IV.9 1 is indeed a mul- 
tiindex. The general method of constructing such an 



operator basis via representation theory is given in Ap- 
pendix VI Of course computing eigenvalues in an ar- 
bitrary basis (also in our basis) is quite hard, because 
complexity of calculations grows very fast with the pa- 
rameter t. 



D. Numerical results 
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0.42486035753 
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5 


0.41022855573 









TABLE I: Numerically calculated spectral gaps for different 
degrees t and qubits number k. 



Here, we present our numerically calculated values 
of spectral gaps for different cases, together with cor- 
responding lower bounds (both '"local"' and '"global"' 
when possible). In Figure |2j there are spectral gaps 
A(Hjtf) for the symmetric group S(2) and different 
number of qubits n, corresponding to the 2-design. 
Both lower bounds are also marked. Similarly, values 
of spectral gaps A(H/ Cf ), for t = 3,4,5 (3, 4, 5-designs) 
can be found in Figures [5J[4] and [5] Notice, that with in- 
creasing t, the '"global"' bound tends to less and less in- 
formation, namely, telling us that gaps are closed. What 
is more, for the last plot (Figure pb, the "local" bound 
does not give any bound since it takes a negative value. 
Knowledge of gaps for k > 4 is required to resolve this 
problem. 

From Figures |2]andj3]one can observe that, in the case 
of t = 2 and t = 3, allspectral gaps are the same. This 
could implies the following 



FIG. 2: (Color online) Numerically evaluated spectral gaps 
Gap = A{H Kt ) for f = 2 (dots) vs. "local" bound for k qubits 
(squares) vs. "global" bound for a given n (here, n = 2) 
(line). 



FIG. 4: (Color online) Numerically evaluated spectral gaps 
Gap = A(Hjy) for t = 4 (dots) vs. "local" bound for k qubits 
(squares) vs. "global" bound for a given n (here, « = 4) 
(line). 



FIG. 3: (Color online) Numerically evaluated spectral gaps 
Gap = A(H kt ) for t = 3 (dots) vs. "local" bound for k qubits 
(squares) vs. "global" bound for a given n (here, n = 4) 
(line). 



Conjecture 8. Knowledge of the spectral gaps of the local 
Hamiltonian from Eq. ( jIH.6) is sufficient to know the 
ones for Hy, since there is a one to one correspondence be- 
tween them. 

Remark 9 Let us explain here one important notation. For 
simplicity and transparency, we mark each lower bound from 
Eq. | |IV.5| corresponding to some spectral gap as a "square" 
in Figures [2] [3] [4] and [5] 



FIG. 5: (Color online) Numerically evaluated spectral gaps 
Gap = A(H^ f ) for f = 5 (dots) vs. "local" bound for k qubits 
(squares). Notice that the '"local"' one does not give any 
bound in this case, because it is negative. Here, the "global"' 
bound can't be reported since a value of a gap for H5 5 is 
needed to calculate it. 



size random quantum circuit is an e-approximate 2- 
design. 

When t = 3, A(H^^) = ^ and thus again the above 
result is valid. The two results are in accordance with 
these from [21 J. 



V. CONCLUSIONS AND OPEN PROBLEMS 



E. Basic examples 



Consider the case when t 
k = 3 is equal to A (^2) = 



2. The spectral g ap for 
Using Theorem 



111.10 



and Fact |5j we have that in this case the second largest 
eigenvalue A2 equals to |. This means that \nlog{\)- 



In this paper, we numerically studied the recent state- 
ment that local random quantum circuits acting on n 
qubits and composed of polynomially many two-qubit 
gates form an approximate unitary poly(n)-design. To 
do this, we evaluated spectral gaps of local Hamilto- 
nians acting on n qubits, using techniques from many- 
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body physics and relating the degree t in f-design to the 
symmetric group S(t). As an additional result, it occurs 
that for t = 2 and 3, all spectral gaps are the same in 
both cases, when there is not any connection between 
the spectral gaps for higher values of t. What is more, 
we compare our results to two lower bounds for spec- 
tral gaps, leading to conclusion that for small number 
of qubits, lower bounding is, usually, not sufficient to 
obtain a reliable result - there is a big difference in ac- 
tual values of spectral gaps and lower bound. For big 
numbers of qubits and high orders of t, there is another 
problem, to obtain a "good" lower bound, one needs to 
calculate first, a spectral gap for a quite big value of n 
and t, which is a quite complicated task, from the com- 
putational point of view. That's why, we were unable 
to compute the "global" lower bound for the case, when 
t — 5, since it requires initial knowledge of the spectral 
gap for k = 5. Here, one possible way, to obtain better 
results would be to use a super-computer and/or the 
power of parallel computing for calculations of spectral 
gaps. An another way, would be to find a "better" ba- 
sis in which our operator X takes a diagonal form or 
at least a block diagonal form. We leave these tasks as 
open. 

At the end, we would like to point out one prob- 
lem which should be of some interest, is it possible 
to approximate values of spectral gaps by some func- 
tion with dependance on parameters k and f? Based on 

the numerics, the function (a(t) 
promising candidate. 

Also, a question, why all spectral gaps for t — 2 and 
f = 3 are the same, remains without an answer. 
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VI. APPENDIX A: ORTHOGONALISATION OF 
REPRESENTATION OPERATORS OF FINITE GROUPS 

In this section we we briefly remind some properties 
of the algebra generated by a given complex finite di- 
mensional representation of the finite group G. The con- 
tent of this section can be found in the standard books 
on representation theory of finite groups and algebras, 
for example, in 13611571 . 

Any complex finite-dimensional representation D : 
G — > Hom(V) of the finite group G, where V is a com- 
plex linear space, generates a algebra Ay[G] C Hom(V) 



r looks like a 



which isomorphic to the group algebra C [G] if the rep- 
resentation D is faithful. Obviously 



A V [G) =span c {D(g), g e G}. 



(VI.l) 



If the operators D(g) are linearly independent, 
then they form a basis of the algebra Ay[G] and 
dim Ay [G] = |G|. It is also possible, using matrix irre- 
ducible representations, to construct a new basis which 
has remarkable properties, very useful in applications 
of representation theory. Below we describe this con- 
struction. 

Notation 10. Let G be a finite group of order \ G\ = n which 
has r classes of conjugated elements. Then G has exactly 
r inequivalent, irreducible representations, in particular G 
has exactly r inequivalent, irreducible matrix representations. 
Let 



D a : G -> Hom(V a ), 



1,2, r, 



dimV a = d a 
(VI.2) 

be all inequivalent, irreducible representations of G and let 
chose these representations to be all unitary (always possible) 
i.e. 



D K (g) = (D?-(g)), 



and (D?.(g)) + = {D«{g))-\ 

(VI.3) 



where i,j = l,2,....,d tt . 

The matrix elements Dfj (g) will play a crucial role in 
the following. 

Definition 11. Let D : G — > Hom(V) be an unitary rep- 
resentation of a finite group G such that the operators D(g), 
g 6 G are linearly independent i.e. dim Ay [G] = |G| and 
let D a : G — >■ Hom(V a ) be all inequivalent, irreducible rep- 
resentations of G described in Notation 10 above. Define 



(VL4) 



1,2,.., d a , Ejj e A V [G] c 



where a = 1,2, ...,r, i,j 
Hom(V). 

The operators have noticeable properties listed in the 

Theorem 12. I) There are exactly \G\ — n nonzero opera- 
tors Efj and 



D(g)=j:DUg)E" 



(VL5) 



l IK 



II) the operators Ef- are orthogonal with respect to the 
Hilbert-Schmidt scalar product in the space Hom{V). 

(E§, Ej) = fr((E«.) + Ef ; ) = k a S^S ik S n , k a > 1, 

(VL6) 

where k a is equal to the multiplicity of the irre- 
ducible representation D a in D and it does not depend on 
i,j = 1,2,...., d u 
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III) the operators Ef- satisfy the following composition rule 

£«4 = 6«%Ef u (VI.7) 

in particular Ef. are orthogonal projections. 

Remark 13 From point II) of the theorem it follows that the 
equations. 



g£G 



(VI.8) 



describe transformation of orthogonalisation of operators 
D(g), g G Gin the space Hom(V) with the Hilbert-Schmidt 

scalar product. 

Remark 14 We can look at operators Efj as a vectors in C"' 
then they are orthonormal (with a coefficient J ^) with re- 



spect to the usual scalar product in C. So in fact we have 
double orthogonality. 

The operators are not only orthogonal projections 
onto their proper subspaces in V but they are also 
orthogonal with respect to the Hilbert-Schmidt scalar 
product in the space Hom(V). 

The basis {Efj} play essential role when D : G — » 

C [G] is the regular representation. In this case the prop- 
erties of the basis {Efj} express the well-known fact that 

the group algebra C [G] is a direct sum of simple matrix 
algebras generated by the irreducible representations of 
the group G. It is always possible to construct the op- 
erators Efj even if the operators D(g) are not linearly 
independent but in this case some of them will be zero. 



From Theorem 12 it follows directly 



Corollary 15. 



r<H12) 



E D ji 

geG 



)D(g)®D(g) 



EEE 

xfi kl mn 



geG 



)Dl k (g)DUg))E%i®rt 



(VI.9) 



thus the coefficients on RHS are expressed only by the matrix 
elements of irreducible representations. 

VII. APPENDIX B: THE LARGEST EIGENVALUE OF 

Gfi,t 

Let us recall one more time the fact that we present 
in Section [ill] 

Fact 16. We have 

\\G F *k it - G^tWco < 4> (vn.i) 

where A2 is the second largest eigenvalue of G^t- Moreover 
the largest eigenvalue Ai of Gu,t is equal to 1, and the corre- 
sponding eigenprojector is equal to Gu H ,t- 

Now, we will prove that the largest eigenvalue of G^t 
is equal to 1. 

Proof. We know that G«f is a operator such that < 
Gu,t _■ !• Having this in mind, let us eigen-decompose 
operator Gu,t in some basis as 

G }t ,t = lE P r,i+l=L V lQl' ( VIL2 ) 

n 1=1 / 

where V\ are eigenvalues, so they follow the constraint 
< V\ < 1, and Q; are the corresponding eigenvectors. 



Now, let us extend Eq. < |VII2) to k walks in the random 
walk: 

since Q/ are projectors. We also have the following 

\mxG tkt = G }lH , t . (VII.4) 

Then, it is easy to observe that the correspondence is 
valid only when V\ — 1 and the eigenprojector Q; is 
equal to G« H/ f. □ 

VIII. APPENDIX C: PROPERTIES OF H, ht 

The Hamiltonian H n r t has the following structure and 
properties Il29ll : 

Lemma 17. For every n, t > the following properties of 
H n ,i = Td(I — Pij+i) hold true: 

1. it is frustration-free: every state \ip) in the groundstate 
manifold, composed of all eigenvectors with eigenvalue 
zero, is such that V/Ff^+i |^) = 0; 

2. the minimum eigenvalue of H Ht t is zero and the zero 
eigenspace is given by 

g n , t :=span{\ip n/d } m , \y nA ) := (Id®V d (7r))|O d ) : n E 

(vni.i) 
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with \<t>j) := d t ^ 1 Y!k = i\k,k) the maximally en- 
tangled state on (C d ) 0t <g> {<C d Y Jt , S t the symmetric 
group of order t, and V d (n) the representation of the 
permutation n G St which acts on (C d )®' as 

v d (7i)\h)®...®\l t ) = 1/^1(1)) ®...®|/„-i( f) >; (vni.2) 



and 



E (^k^^kA) 

neSt 



< 



2? 
d" ' 



(VIII.4) 



3. let G nr t be the projector onto Q„ )t . Then 

2t 2 

E i<vvii/vd)r<i + ^ 

TTGSf " 



(VIII.3) 
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